Thermodynamics of a gas of deconfined bosonic spinons in two dimensions 
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We consider the quantum phase transition between a Neel antiferromagnet and a valence-bond 
solid (VBS) in a two-dimensional system of S = 1/2 spins. Assuming that the excitations of the 
critical ground state are linearly dispersing deconfined spinons obeying Bose statistics, we derive 
expressions for the specific heat and the magnetic susceptibility at low temperature T. Comparing 
with quantum Monte Carlo results for the J-Q model, which is a candidate for a deconfined Neel- 
VBS transition, we find excellent agreement, including a previously noted logarithmic correction 
in the susceptibility. In our treatment, this is a direct consequence of a confinement length scale 
A <x £ 1+a oc l/T 1+a , where £ is the correlation length and a > (with a ~ 0.2 in the model). 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Mg, 75.40.Cx 



The excitations of a quantum antiferromagnet are spin 
waves (magnons) carrying spin S = 1. At a conventional 
quantum phase transition in which the antiferromagnetic 
(Neel) order parameter vanishes continuously 0, Q , the 
magnons remain well-dehned elementary excitations even 
at the critical point. The gapped "triplon" excitations in 
the nonmagnetic phase also have S — 1. A different sce- 
nario has also been proposed, in which the magnons of 
a two-dimensional system fractionalize into independent 
(deconfined) S = 1/2 spinons at the phase transition, 
which was termed the deconfined quantum critical (DQC) 
point Q. In the non- magnetic phase, which in this case 
is a valence-bond solid (VBS) with spontaneously bro- 
ken lattice symmetries [J], the spinons are confined into 
excitations which carry spin S — 1 and 5 = 0. This 
proposal is supported by quantum Monte Carlo (QMC) 
simulations of a "J-Q" model [5j-|9| (an S =1/2 Heisen- 
berg model including four-spin or higher-order terms), 
but objections have also been raised [10h12|. The decon- 



finement scenario violates the long-held "Landau rule" 
according to which order-order transitions breaking un- 
related symmetries should be of first order. 

In this Letter we provide evidence for deconfinement in 
the J-Q model based on its thermodynamic properties. 
Using a phenomenological ansatz of deconfined bosonic 
spinons with linear dispersion e(fc) = ck at temperature 
T = 0, we derive expressions for the specific heat and 
the magnetic susceptibility at T > 0. In addition to the 
spinon velocity c, these quantities depend on the length 
scale A(T) within which spinons are confined at T > 0. 
As a consequence of the T — DQC theory [§] , one may 
expect this length to diverge when T — > as a power of 
the correlation length, A cx £ 1+a , where the anomalous 
exponent a > and the correlation length diverges as 
£ oc 1/T Z with dynamic exponent z = 1. The correlation 
length is defined in the standard way, using the antiferro- 
magnetic or VBS correlations (which diverge in the same 
way at the DQC point). We will show that confinement 



on the larger length scale when a > implies a logarith- 
mic correction to the standard quantum-critical scaling 
X oc T of the susceptibility, which agrees with a recent 
QMC study of the J-Q model Q- The specific heat and 
the susceptibility obey the spinon-gas forms with a com- 
mon value of c and a « 0.2. 

Spinon gas — Our assumption is that a system with 
couplings tuned to the T = quantum-critical point can 
be described as a gas of bosonic spinons with dispersion 
e(k) = ^Jc 2 k 2 + A 2 (T) at T > 0. This dispersion is valid 
for magnons at a conventional 0(3) quantum phase tran- 
sition between the Neel state and a disordered state (e.g., 
in dimerized Heisenberg models fl3l - fl9j |). in which case 
the thermal "gap" A is related to the correlation length 
£ according to A oc l/£ oc T z , with z = 1 Q. At 
the DQC point £ oc 1/T is also expected, but the spinon 
gap should be given by the larger confinement length, 
A oc 1/A oc T 1+a . Effectively, the gap is used as an 
infrared cut-off in momentum space. 

The T = confinement exponent a was previously 
estimated using QMC results for the finite-size scaling 
of the U(l)-Z4 cross-over of the VBS order-parameter 
symmetry (a hallmark of the DQC theory [3]) of a variant 
of the J-Q model. The result was a = 0.20 ± 0.05. It is 
unclear whether the T > confinement exponent should 
be the same, however. Here we leave a as a free parameter 
along with the velocity c and write the gap as 



mi/ 2 T (T/c) a . 



(1) 



where the constant mi/2 should be close to 1. In the 
case of magnons, it is known that A = m,\T ', with the 
mean-field value mi ~ 0.96 [2| in good agreement with 
QMC calculations of observables (e.g., the magnetic sus- 
ceptibility) which depend on this constant 2. 161. 
In a magnetic field B, the spinon level is split into 



e± (fc) = \/c 2 k 2 + A 2 ± nB = e{k) ± fj,B, (2) 
where fj, = 1/2. This form with fi = 1 holds also for 
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the two shifted magnon levels (with eo not shifted). In 
the CP 1 DQC theory Q, there are both spinons and 
anti-spinons, which contribute equally to thermodynamic 
properties. We take this into account with a factor F = 2, 
while for magnons F = 1. 

With the boson occupation number n(e) = l/(e c / T — 1) 
the magnetization per lattice site for small B is: 



M = fiF 
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lo sinh 2 [i^/a ; 2 + (A/T) 2 ] 
The integral can be computed exactly, 
xdx 4p 



(3) 



smh 2 (^x 2 + p 2 ) 



1 - e~P 



41n(e p -l), (4) 



where p = A/T. For magnons at the usual 0(3) transi- 
tion, p = mi « 0.96 and the susceptibility is [2J 



Xi « (1.0760/ttc 2 )T. 



(5) 



For spinons, if there is indeed an anomalous exponent 
a > in Eq. (JTJ, then A/T ->OasT^0 and we can 
use the expansion of ((3]) around p = 0, giving 
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(6) 



Here we have used m\/2 = 1 in ([T]), and the next cor- 
rection to X1/2/T is of order (T/c) 4a . The logarithmic 
correction is very interesting, as it was already identified 
in a recent QMC study of the J-Q model 0. 
The specific heat per site is 



C s = (2S + 1)F / e(k) 



dn{e) d 2 k 



dT (2tt) 2 ' 

which for S = 1/2 leads to the low-T behavior 



2T z 
C1/2 = T X 



7TCT 



6C(3) 



(7) 



(8) 



+ a + a(l + a) hi 



where we have again used mi/2 = 1 and £(3) ~ 1.20206. 
Note that the log-correction is multiplied by a power and 
is not as dramatic as in the susceptibility ©. For the 
0(3) transition, S = F = 1 in © gives, at low T p|.[20|. 



Ci = [36C(3)/5ttc 2 ]T 2 



(9) 



Apart from the logarithms, the differences in the ther- 
modynamics between spinons and magnons arise mainly 
from the degeneracy factors and fi. The log correction 
to X1/2 in © is significant, however, and should be a 
decisive fingerprint of the deconfined spinon gas. 
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FIG. 1: (Color online) Size dependent critical coupling of the 
J-Q model based on (L, 2L) curve crossings of the scaled spin 
and dimer correlation lengths, £ S /L and £d/L. The average of 
these two estimates is also shown. The horizontal line shows 
the critical point extracted from the average. 

Quantum- critical models — A promising model exhibit- 
ing a Neel-VBS transition is the J-Q model, which in its 
simplest form is defined by the hamiltonian [j| 



H 



(ijkl) 



kl, 



(10) 



where is a singlet projector; CV,- = 1/4 — Sj • Sj. 
In the J (Heisenberg) term ij are nearest neighbors on 
the square lattice, while in the Q term ij and kl form 
opposite edges of a 2 x 2 plaquette. There is mounting 
evidence [5|-|9j of a continuous T = transition in this 
system between a Neel state for J/Q > (J/Q) c an d a 
VBS for J/Q < (J/Q) c , with (J/Q) c 0.045 Q. 

To ensure pure critical behavior when testing the above 
spinon gas predictions, it is necessary to know the criti- 
cal coupling ratio ( J/Q) c to high precision. Here we use 
the correlation lengths £ s and £<j extracted from, respec- 
tively, the spin-spin and dimer-dimer (four-spin) corre- 
lation functions. QMC calculations were carried out on 
L x L lattices with L up to 256 at inverse temperature 
P = Q/T = L (as in Q). At a DQC point, £ S /L and 
£,d/L should be size independent for large L. Curves 
plotted versus J/Q for two different system sizes, e.g., 
L and 2L, should then cross each other at some value 
(J/Q)l, which can be different for £ S /L and £<j/£ but 
in both cases should approach {J/Q) c when L — > 00. 
Such crossing points are shown in Fig. [TJ Due to the 
slow convergence, it is difficult to extrapolate precisely. 
There is, however, a remarkable feature of these data: 
The £ S /L and S,d/L crossing points approach an appar- 
ent common asymptotic value at the same rate but from 
different sides. Their average exhibits almost no size 
dependence, and one can therefore obtain a much bet- 
ter critical-point estimate than what might initially have 
been expected. The result based on the four largest-L 
points (which agree completely within statistical errors) 
is {J/Q)c = 0.04498(3). Here we will use J/Q = 0.045. 



3 



0.28 



0.26 



til 0.24 
ESS 



0.22 



0.20 



J-Q model 




FIG. 2: (Color online) Temperature dependence of the energy 
E relative to the ground state energy Eq with the expected 
leading T 3 dependence divided out. The horizontal lines show 
the prefactor of the cubic term in the E(T) fits. 



Conventional 0(3) T > scaling has been studied in 
the past in various dimerized Heisenberg models (where 
the hamiltonian itself breaks lattice symmetries, and no 
other symmetries are broken in the disordered phase) 
l44l7||. To compare with the J-Q model, we consider 
a system with couplings J and J' > J, with the stronger 
ones arranged in columns. This model was the subject 
of a recent high-precision study [19j , which gave the crit- 
ical ratio (J'/J) c = 1.9096(2). A further improved esti- 
mate is now available, (J'/J) c = 1.90948(4) We use 
J' / J = 1.9095 and carry out more detailed comparisons 
with the 0(3) theory than in past studies. In this case 
all quantities should be normalized per two-site unit cell. 

Energy fits — We first fit low-T results for the internal 
energy based on the leading specific heat forms ((5J and 
© . QMC calculations were carried out using sufficiently 
large lattices (L < 512) to eliminate finite-size effects in 
the range of temperatures considered. We also use the 
ground state energy extrapolated to L = oo based on 
(3 = L results. For the J-Q model at J/Q = 0.045, the 
result is E a /Q = -0.8740318(4), while the J- J' model 
with J'/J = 1.9095 has E /J = -1.740507(2) per unit 
cell. Fig. [5] shows E(T) after Eq has been subtracted 
and T 3 has been divided out. The low-T behavior gives 
the spinon velocity c = 2.55Q for the J-Q model and 
the magnon velocity c = 1.38J for the J-J' model (which 



should be interpreted as c = y/c x c y since the J-J' model 
is anisotropic). While there are corrections to the T 3 
behavior in Fig. [5J the low-T results for the J-Q model 
are not sufficiently accurate to test the correction in © . 
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FIG. 3: (Color online) Susceptibility divided by T. The 
dashed curve in the J-Q graph is Eq. ((6| with c = 2.55Q, 
a — 0.22 (without the power-law correction). The solid curve 
has c = 2.48Q, a = 0.20, and also includes a T 2 term. The 
solid curve in the J-J' graph is a second-order polynomial. 
The insets show the same data on logarithmic T scales. 

It is anyway doubtful whether the spinon gas model can 
correctly capture subleading corrections. 

Susceptibility fits — The susceptibilities of both models 
are shown in Fig. [3] with T divided out. In the case of 
the J-J' model, there are significant corrections to the 
asymptotic T — > constant behavior expected with (JS). 
A second-order polynomial fit to the low-T data is shown. 
The extrapolated T = susceptibility corresponds to a 
spin- wave velocity c = 1.36 J in ([5]), in excellent agree- 
ment with the value obtained from the energy. 

Two different fits to J-Q results are shown in Fig. [3] 
Neglecting the power-law correction in ((5]) and fixing the 
spinon velocity to the value c = 2.55Q obtained from the 
energy, only the exponent a in front of the log correction 
was adjusted to match the low-T data. This result is 
a = 0.22. A completely independent fit to a wider range 
of data, and including a T 2 correction (which was found 
to describe the data well in [8j), gives c = 2A8Q and 
a = 0.20. In both cases, including also the very small 
power-law correction in ((SJ changes X1/2/T by less than 
1% and barely affects the extracted parameters. 

We have set mx/2 = 1 in © throughout the above 
analysis, while we may only expect mi/2 ~ 1 [as with 
the constant mi in the conventional 0(3) theory]. Phys- 
ical observables depend only weakly on m-1/2, however, 
and the consistent c-values extracted from two different 
quantities justify the use of — la posteriori. 

Wilson ratio — The Wilson ratio of the J-Q model ex- 
hibits a log divergence. From Eqs. © and © and the 
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parameters of the fits (taking c from the energy fit and 
a from the susceptibility fit with c fixed) we get W1/2 = 
XT/C = w 1/2 [l+ala(c/T)}, with w 1/2 = 0.0346±0.0002, 
c = 2.55±0.02, and a = 0.222±0.005. For the J-J' model 
we get Wi — 0.f262 ± 0.0006, in good agreement with 
Wi = 0.1243 from Eqs. © and ©. Including the next 
term in the 1/N expansion of the large- N 0(3) theory 
[3 makes this agreement worse by several percent, how- 
ever. Note that if the log-correction is disregarded, Wi/ 2 
is only about 1 /4 of W\ . 

Conclusions and discussion — We have tested a model 
of non-interacting (deconfined) bosonic spinons against 
QMC data for the J-Q model, which is a promising can- 
didate for a DQC point. The most notable result is that 
a confinement length A diverging as l/T 1+a with a > 
leads to a logarithmic correction to the susceptibility \i 
as was previously observed in the J-Q model [8J. The ve- 
locity entering in \ agrees with the velocity needed to fit 
the specific heat. The anomalous exponent a s» 0.22 is in 
good agreement with an estimate based on a completely 
different analysis at T = Q, which suggests that the 
T = and T > exponents indeed are the same (which 
is unclear in the DQC theory, in which no anomalous 
T > exponent has been discussed d S3). The criti- 
cal behavior does not fit the standard 0(3) picture with 
5 = 1 excitations (3, which we have investigated here in 
the context of a dimerized model. 

The agreement between the critical J-Q model and 
the non-interacting spinon gas is remarkable, considering 
that the spinons in the DQC theory are only marginally 
deconfined (with interactions mediated by the gauge 
field) (3- Apparently, beyond their underlying role in 
determining the anomalous exponent a in Eq. |T]), these 
interactions only have very small effects on the thermo- 
dynamics. A treatment similar to the spinon gas consid- 
ered here has been applied to the 5=1/2 Heisenberg 
chain (with the important difference that the spinons 
there obey Fermi statistics) [23|. Known results, includ- 
ing logarithmic corrections, were reproduced. 

It would be useful to have an independent estimate of 
the spinon velocity. A velocity c = 2.4±0.3 was extracted 
for the critical J-Q model in Q, using a criterion for cu- 
bic space-time geometry in QMC simulations. Although 
the value is in good agreement with ours, it is unclear 
whether their method applies to spinons (while it should 
work for magnons). In future studies we will extract c 
from imaginary-time dependent spin-spin correlations. 

Our study lends support to the DQC scenario [3 for 
the Neel-VBS transition, although the phenomenological 
approach does not address the mechanism of deconfine- 
ment (only tests the consequences). Log corrections at 
T > should also have counterparts at T = 0. Further 
work along these lines will hopefully explain, e.g., anoma- 
lous corrections to the spin stiffness of the J-Q model [3 
and its impurity response [9( . An important missing link 
is how these corrections could arise from the CP 1 field 



theory of the DQC proposal [3, i.e., whether this theory 
is complete in its current form or whether some ingre- 
dient is still missing. No log corrections were found in 
large- TV treatments of the CP^" 1 theory [I2,[24|], but it 
is possible that these corrections appear only for small 
N. A logarithmic enhancement of the susceptibility was 
found in a U(l) gauge theory with fermions [25(. In that 
case, there is also a correction to the specific heat, which 
makes the Wilson ratio non-divergent. The spinon gas 
approach with Fermi statistics gives no log corrections 
and the Wilson ratio equals 0.0320, which is identical 
(with fi = 1/2) to the value in Ref. [22 . 



We finally note that there are no indications of a first- 
order transition in the J-Q model (with previous claims 
fiol 11 1 not supported by later results [1, (3)- As a 
matter of principle, however, extremely weak discontinu- 
ities cannot be ruled out based on numerical data alone 
(though the first-order scenario appears increasingly un- 
likely). What we have shown here is that, regardless of 
the ultimate nature of the transition, spinons are decon- 
fined on length scales sufficiently large to have significant 
consequences for the thermodynamics. 
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